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Abstract 



We present a general result within the bidifferential calculus approach to integrable partial 
differential and difference equations, which can be regarded as a quite universal formulation of 
a vectorial binary Darboux transformation. This is then applied to the non-autonomous chiral 
model, a certain reduction of which is known to appear in the case of the D-dimensional vacuum 
Einstein equations with D — 2 commuting Killing vector fields. A large class of exact solutions 
■ is obtained, and the aforementioned reduction is implemented. This results in an alternative to 

1 the well-known Belinski-Zakharov formalism. We recover relevant examples of space-times in 

dimensions four (Kerr-NUT, Tomimatsu-Sato) and five (single and double Myers-Perry, black 
saturn, bicycling black rings). 

; 1 Introduction 

O 

The bidifferential calculus formalism aims to understand integrability features and solution gen- 
erating methods for (at least a large class) of integrable partial differential or difference equations 
(PDDEs) resolved from the particularities of examples, i.e. on an as far as possible universal level 
[1-3] . The most basic ingredient is a graded associative algebra, supplied with two anti-commuting 
graded derivations of degree one. It can and should be regarded as a generalization (in the spirit of 
noncommutative geometry) of the algebra of differential forms on a manifold, but supplied with two 
analogs of the exterior derivative. Once a PDDE is translated to this framework, it is extremely 
simply to elaborate its integrability conditions. In fact, it could not be simpler. 

In this framework, Darboux transformations (see [4-6] and the references therein) have been first 
addressed in [3]. In the latter work, we had also obtained a very simple solution generating result 
that evolved into improved versions in more recent applications [7-9]. Its relation with Darboux 
transformations has been further clarified in [9] (see Appendix A therein), an essential step toward 
the much more general result that we present in this work. The resulting class of solutions is 
expressed in a universal way, in the sense that the corresponding formula holds simultaneously 
for all integrable PDDEs possessing a bidifferential calculus formulation. Choosing a bidifferential 
calculus associated with a specific PDDE, we can generate infinite families of (soliton-like) solutions. 

In Section 2 we recall some basics of bidifferential calculus and set up the stage in a concise 
way to formulate our general result about Darboux transformations. This hardly requires previous 
knowledge and can be taken as an independent and self-contained step into the world of integrable 
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PDDEs (though perhaps leaving aside a lot of other fascinating aspects of the latter). That this is 
indeed a very powerful tool, is demonstrated in the rest of this work, where we concentrate on one 
of the more tricky examples of integrable PDEs. 

In Section 3 we elaborate in detail the example of the non-autonomous chiral model equation 

(P9z9~ 1 )z + e(pg p g- 1 ) p = (1.1) 

for an m x m matrix g, where p > and z are independent real variables and e = il. 1 With e = 1, 
this governs the case of the stationary, axially symmetric vacuum Einstein (m = 2) and Einstein- 
Maxwell (m = 3) equations in four dimensions, where we have two commuting Killing vector 
fields, one spacelike and the other one asymptotically timelike (see e.g. [10-13]). In an analogous 
way, (1.1) with m > 3 appears in the dimensional reduction of (the bosonic part of) higher- 
dimensional supergravity theories to two dimensions (see e.g. [14-18]). With e = —1, the above 
equation appears in the case of vacuum solutions of Einstein's equations in four dimensions with 
two commuting spacelike Killing vector fields [11,19-21], describing e.g. cylindrical gravitational 
waves. 

In Section 3, we also present a reduction condition that, imposed on the obtained family of 
solutions, achieves that g is symmetric. Any such real and symmetric g determines a solution of 
the vacuum Einstein equations inm + 2 dimensions, as recalled in Section 4. The resulting recipe 
to construct solutions of the vacuum Einstein equations is close to the familiar Belinski-Zakharov 
method [10,11,19], which has been applied in numerous publications. In a sense, what we obtained 
is a kind of matrix version of the latter. 

In Section 4 we work out several examples in four and five space-time dimensions, and show 
that the resulting metrics include important solutions of Einstein's equations, in particular some 
of the more recently found black objects in five dimensions, which have no counterpart in four 
dimensions (see e.g. [16,18,22]). Here we used Mathematica [23] in a substantial way. 

Section 5 contains some concluding remarks. 



2 Darboux transformations in bidifferential calculus 

A graded associative algebra is an associative algebra CI over C with a direct sum decomposition 
ft = r>o ^ r into a subalgebra A := ^° and .A-bimodules Cl r , such that Cl r ft s C ft r+s . A 
bidifferential calculus (or bidifferential graded algebra) is a unital graded associative algebra f2, 
equipped with two (C-linear) graded derivations 2 d,d : ft — > $7 of degree one (hence dSl r C $T +1 , 
d$T C $T +1 ), with the properties 

d 2 = d 2 = dd + dd = 0. (2.1) 

In the bidifferential calculus approach to integrable PDDEs we are looking for a choice (fi, d, d) 
and some <f> € A, or some invertible g G A, such that either 

dd0 = d</>d0 (2-2) 

or 

d[(dg)g- 1 }=0 (2.3) 

X A subscript indicates a partial derivative with respect to the corresponding variable. Formally, p M> ip relates 
the two values of e in (1.1). 

2 Hence d and d both satisfy the graded Leibniz rule d(xx') — (dx) X + ( — l) r X&X > f° r au X £ ^ r an d X £ £1- 
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is equivalent to a certain PDDE. The two equations are related by the Miura equation 

(dgjg-^dtp, (2.4) 

which establishes a kind of Miura transformation between the two equations (2.2) and (2.3). This 
equation has both, (2.2) and (2.3), as integrability conditions. This means that if we find a solution 
pair (<j),g) of (2.4), then <fi solves (2.2) and g solves (2.3). 

In the following theorem, we formulate a solution-generating result which amounts to a trans- 
formation that takes a given solution (</>o,<7o) of the Miura equation (2.4) to a new solution of 
(2-4), 

{4>o,9o) ^ (fag) ■ 

This induces corresponding transformations of solutions of (2.2), respectively (2.3). 

Let now A be the algebra of all finite-dimensional 3 matrices, with entries in a unital algebra B. 
The product of two matrices is defined to be zero if the sizes of the two matrices do not match. In 
the following we assume that there is a graded algebra Q. with Q° = A, and a bidifferential calculus 
($7,d,d), and such that d and d preserve the size of matrices. I = I m and I = I n denote the 
m x m, respectively n x n, identity matrix, and we assume that they are annihilated by d and d. 
Furthermore, Mat(m, n, £>) denotes the set ofmxn matrices over B. 

Theorem 2.1. Let (fio,go G Mat(ra, m, £>) solve the Miura equation (2.4)- Let P,Q G Mat(n,n,S) 
be invertible solutions of 

dP = (dP)P, dQ = QdQ, (2.5) 

which are independent in the sense that QY = Y P implies Y = 0. Let U G Mat(m, n,B) and 
V G Mat(n, m, B) be solutions of the linear equations 4 " 

~dU = (dU) P + (d^o) U , dV = QdV- V d<p ■ (2.6) 

Furthermore, let X £ Mat(n,n, B) be an invertible solution of the Sylvester-type equation 

X P - QX = VU . (2.7) 

Then 

<j ) = ^ + UX- l V, g = (L + U(QX)- 1 V)g (2.8) 
solve the Miura equation (2.4), and thus also (2.2), respectively (2.3). 

Proof. Acting on (2.7) with d, using the rules of bidifferential calculus, (2.5), (2.6), and (2.7) again, 
we find 

Q [dX - (dX) P + (dQ) X + (dV) U] = [dX - (dX) P + (dQ) X + (dV) U] P . 
Since P and Q are assumed to be independent, this implies 5 

dX-(dX)P+(dQ)X + (dV)U = 0. (2.9) 

3 An extension to a suitable class of oo-dimensional matrices, respectively operators, is certainly possible. 

4 If we are primarily interested in solving (2.3), it is more convenient to replace dcj>o by (dgo) g^ 1 in these equations, 
by use of the Miura equation for (0o,<?o)- 

5 It is a key feature of Theorem 2.1 that this equation is automatically solved if P and Q are independent. Relaxing 
the latter assumption, the theorem remains valid if this equation is added to the assumptions. 
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Using also (2.5) and (2.7), we easily deduce that 

diQxy 1 = -X- 1 [(dX)X- 1 (XP)- (dVWiQX)- 1 

= (dX- 1 ) [I + VU (QX)- 1 ] + X-^dV) U(QX)- 1 . 

This is then used in the elaboration of 

dg = d[(I + U(QXr 1 V)g } 

via the graded derivation (Leibniz) rule for d. We eliminate all further terms involving a d with 
the help of (2.6), and apply (2.7) to eliminate a P in favor of a Q. Finally we obtain dg = (d<p) g 
with cfr given by (2.8). Our assumptions ensure that the inverse of g exists. It is given by 

g- 1 =9o\l-U{XP)- x V). 

The integrability conditions of (2.6) and (2.9) are satisfied as a consequence of (2.5) and dd0 o = 
d</>od0o (which follows from (dg$) g^ 1 = dcfro). □ 

The similarity with corresponding results about (vectorial) binary Darboux transformations in 
the literature on integrable systems is evident (see e.g. [4-6]). As already mentioned in the intro- 
duction, Theorem 2.1 considerably generalizes previous, more restricted results, see in particular 
[9]. 

Remark 2.2. As a consequence of (2.5), P and — Q solve the n x n version of (2.2). Since we 
assume that they are invertible, P and Q 1 solve the n x n version of (2.3). If U is a solution of 
the first of (2.6), then also UP, by use of the first of (2.5). If V is a solution of the second of (2.6), 
then also QV, by use of the second of (2.5). 

Remark 2.3. By direct computations one verifies the following "scaling property". Let U and V 
solve (2.6) and let 

g = wg , U = UW 1 , V = W 2 V, 

with a solution (f,w) of the m = 1 version of (2.4), and invertible Wj 6 Mat(n, n, B), subject to 

[W 1 ,P} = 0, [W 2 ,Q} = 0. 

We further assume that ip and w are in the center of O. 6 If W\ and W 2 satisfy 

dWi = (dWi) P + (dw) w' 1 Wi , dW 2 = Q dW 2 - W 2 (dw) ur 1 , 

then U and V solve (2.6), with </>o replaced by cpo = (fro + ip I. Furthermore, X = W 2 XW\ solves 
(2.7) with U,V replaced by U, V. For the new solutions of (2.4) (and thus (2.2), respectively 
(2.3)), determined by the theorem, we find 

4> = 4> + (pI , g = wg, 

where <fr and g are given by (2.8). Hence, multiplication of g$ by a "scalar" solution w of (2.3) 
simplify amounts to multiplication of g by w. In the special case where w is the identity ele- 
ment and g thus remains unchanged, the above equations for W\ and W 2 still allow non-trivial 
transformations of U and V. These results will be used in Section 4. 

6 This assumption will be a (perhaps too) severe restriction in case B involves differential of difference operators. 
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Remark 2.4. Let (4>i,9i), i = 1)2, be two solutions obtained via Theorem 2.1 from the same seed 
solution (</>o,<7o)- Let (Pj, Qj, [/j, Vj) be corresponding solutions of (2.5) and (2.6). Then 

also solve (2.5) and (2.6). If -P and Q are independent and if (2.7) has an invertible solution X, 
then (2.8) determines a new solution of (2.4), and thus new solutions of (2.2) and (2.3). This 
expresses a nonlinear superposition principle. By iteration, one can build superpositions of an 
arbitrary number of "elementary" solutions, hence (analogs of) "multi-solitons" . 

Remark 2.5. We note that the first of equations (2.6) is a special case (and n x n matrix version) 
of the general linear equation 

d^-A^ = d^P where A = respectively 

(dfiO 9 ■ 

Here P is a solution of the nonlinear equation 

dP = (dP) P . 

An n x n matrix version of it appears in (2.5). The integrability condition of (2.10) is (2.2), 
respectively (2.3). A similar statement holds for the second of equations (2.6), the "adjoint linear 
system" . 

Remark 2.6. In the special case where O is the algebra of matrix- valued differential forms on a 
manifold, and d the exterior derivative, any tensor field M of type (1, 1) with vanishing Nijenhuis 
torsion determines a map d_^ satisfying the above conditions [24, 25]. Moreover, according to 
Frolicher-Nijenhuis theory [24], any d such that (2.1) holds has to be of this form. Finite-dimensional 
integrable systems have been considered in this framework in [26]. The generalization to Lie 
algebroid structures is nicely described in [27]. See also [28-30] for related aspects. We should 
stress, however, that in our central examples we depart from differential geometry and consider a 
differential calculus in a weaker sense (e.g. of noncommutative geometry). 

The setting of the above theorem allows in principle structures far away from classical cal- 
culus and in particular differential geometry, then dealing with equations beyond differential and 
difference equations. A somewhat more restricted framework is given by setting 

n = A®/\(c N ), 

where f\(C N ) denotes the exterior (Grassmann) algebra of the vector space C^. In this case it 
is sufficient to define suitable operators d and d on A, since they extend to 0, in an evident way. 
Many integrable PDDEs have been treated in this framework and in the next section we turn to 
an important example. 

3 Solutions of the non-autonomous chiral model 

We consider the bidifferential calculus determined by 

df = -f z 6 + e e (f p - p-'fe) 6 , df = e- e (f p + p~ 1 f e ) fr + ef, & (3.1) 
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for / G C°°(R 3 ) (cf. [9]). 7 £i,£ 2 is a basis of ^(C 2 ). d and d extend via d(/i£i + / 2 £>) = 
(d/i) A 6 + (d/ 2 ) A £2 and d(/& A 6) = (d/) A & A 6 = to .A® A(C 2 ), with B = C°°(M 3 ). For 
an m x m matrix-valued function 5, (2.3) now takes the form 

{P9z9~ x )z + t{pg P g~ 1 )p - [(g P + p~ 1 ge)g~ 1 \e + (geg~ 1 ) P = o, 

which reduces to (1.1) if g does not depend on 9. The coordinate 9 is needed to have the properties 
of a bidifferential calculus, but addressing (1.1) we are primarily interested in equations for objects 
that do not depend on it. 

Using this bidifferential calculus, the Miura equation (2.4) decomposes into 

{g P + P^ge) g' 1 = -e e 4> z , g z g' 1 = e e e (<p p - p~ x <\>e) ■ 
If g is ^-independent, this requires 

<j) = e- 4>, (3-2) 
with ^-independent (f>, and then reduces to 

g P g~ 1 = -(f>z, tg z g~ x = 4> P + p~ x <\> • (3.3) 

In the following, we elaborate Theorem 2.1 using the above bidifferential calculus. Section 3.1 
provides the complete solution of the two nonlinear equations (2.5) under the condition that P 
(respectively Q) has geometrically simple spectrum (i.e. for each eigenvalue there is a unique 
corresponding Jordan block in the Jordan normal form), also see [9]. In this case the two equations 
actually coincide. Then it only remains to solve linear equations, see Section 3.2 and Section 3.3 
below. In Section 3.4 we present a condition to be imposed on the data in order to achieve that 
the solution g is symmetric (or Hermitian). This reduction is crucial in the context of Einstein's 
equations, see Section 4. Most of the following is, however, independent of additional assumptions 
and provides a general procedure to construct solutions of the non-autonomous chiral model. The 
corresponding equations cannot be solved explicitly without some restrictions, in particular on the 
form of the seed solution. We content ourselves with providing illustrative and important examples 
relevant in the context of gravity in Section 4. 

3.1 The equations for P and Q 

In terms of 

P := e e P, 

the first of the two equations (2.5) decomposes into the following pair of equations, 

P p + p-\P e -P) = -P Z P, eP z = [P p - p-\P e -P)]P, (3.4) 

~ 2 

which are autonomous in the variable 9. Assuming that P and I + eP are invertible, and that P 
is ^-independent, (3.4) implies 

p 2 -2p~ l {zl- A)P-eI = 0, (3.5) 

with an arbitrary constant n x n matrix A (also see [9]). This is a matrix version of the pole 
trajectories in the Belinski-Zakharov approach [11,19]. A well-known symmetry of the latter extends 

to the matrix case: (3.5) is invariant under P H> —eP . For the following result, see Lemma 4.1 
m [9]- 

7 The non-autonomous chiral model is well-known to be a reduction of the (anti-) self-dual Yang-Mills equations, 
for which a very simple bidifferential calculus exists [3] (that may be considered as a prototype). From the latter one 
can then derive the bidifferential calculus under consideration. 
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Lemma 3.1. Any 6 -independent, invertible solution of (3.5), which commutes with its derivatives 
with respect to p and z, solves (3.4)- O 

A can be taken in Jordan normal form A = block-diag(A ni , . . . ,A 1ls ), without restriction of 
generality, and a solution of (3.5) is then given by 

P = block-diag(P ni , . . . , P n J , 

where the block P ni is a solution of (3.5) with A replaced by the Jordan block A ni , see the next 
examples. Under the assumption that P has geometrically simple spectrum, this is the most general 
solution of dP = (dP) P with ^-independent P = e 8 P [9]. 

Example 3.2. If A is diagonal, i.e. A = diag(ai, . . . , a n ) with constants a«, corresponding solutions 
of (3.5) are given by P = diag(pi, . . . ,p n ), where pi is any of 



Pi = p 1 (z - a,i + %,) , pi = p 1 (z-a.i-%), 



(3.6) 



with 9\i = \J {z — ai) 2 + e p 2 . Note that pi = —e/p-i- 

Notation 3.3. For better comparison with the relevant literature, we will sometimes write 



P 

Pi = — 

m 



p 

Pi = — , 
Pi 



(3.7) 



where 



Pi 



V(-z - ai) 2 + ep 2 



[z - ai 



Pi 



V (z - ai) 2 + ep 2 - (z - ai) . 



(3. 



Pi and pi are often referred to as (indicating the presence of) a soliton and an anti-soliton, respec- 
tively. Note that pi = — e p 2 / pi. If e = 1, then pi is non- negative and only vanishes on a subset of 
{P = o}. 



Example 3.4. For an r x r Jordan block 



^4.^ Qj 1 y ~\~ ./\^ f , 



(3.5) has the solutions (see [9]) 



N r = 



( 1 
1 



V" 



1 



Pr = p- 1 (zlr-Ar + Y, (i^) 1 ^ [2(a - z) N r + A^) , 



where *R = \/{z — a) 2 + ep 2 . This is an upper-triangular Toeplitz matrix and thus commutes with 
its derivatives. In particular, we have P± = p = p~ x [z — a ± £H] and 



2 = V 



1 =F^H 
1 



-i 



1 t^ _1 ±§/9$K~ 3 p _1 
P 3 =p | 1 1 
1 
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For the solutions P (and thus P), obtained in this way, we have 

dP = (dP) P = P (dP) , 

so that they also provide us with solutions Q = e~ 9 Q, with ^-independent Q, of the second of 
(2.5), also see [9]. 

Recall that on the way to the above results we assumed that P and I + eP are invertible, and 

~ 2 

then also Q and I + eQ . These assumption will also be made throughout in the following. 

Remark 3.5. The source matrix A for P and the corresponding source matrix A' for Q can be 
assumed to be simultaneously in Jordan normal form, without restriction of generality. Since this 
is achieved by similarity transformations with constant transformation matrices, the latter can be 
absorbed by redefinitions that restore all the equations obtained from Theorem 2.1. But in general 
this will no longer be so if we impose a reduction condition that relates P and Q, as in Section 3.4 
below. 

3.2 The equations for U and V 

Setting 0o = e~ 6 (f>o with a ^-independent <fto (cf. (3.2)), (2.6) is autonomous in the variable 6. 
Assuming that U, V are ^-independent, and using the fact that 4>q has to solve the Miura equations 
(3.3) together with some go, we obtain the following systems of linear differential equations for U 
and V, 

U p = (go, p g - 1 U-g , z g G - 1 UP)(I + eP 2 )- 1 , 

U z = (g , zgo - l U + eg , p g l UP) (I + eP 2 )- 1 , (3.9) 

and 

V p = (I + eQY^QVgo^g^-Vgo^gu 1 ), 

V z = -(I + eQ^iVg^g^ + eQVg^g^). (3.10) 

These equations have to be solved for the given seed solution go of (2.3). For diagonal go, this is done 
in the next example for the ^-equations. Similar results are easily obtained for the [/-equations. 

Example 3.6. Let 

Q = diag(pi, . . . ,p n ) , g = diag(wi, . . . , w m ) , 

where pi is either pi or pi in (3.6), and w a is a non- vanishing solution of the scalar (i.e., m = 1) 
version 

[p (In w)p] p = -e [p(\nw) z ] z (3.11) 

of the non-autonomous chiral model (1.1). Writing V = (Via), where i = 1, . . . , n and a = 1, . . . , m, 
(3.10) reads 

(\nV ia ) p = 1 _ 2 (pi (\nw a ) z - (\nw a ) p ) , (\nV ia ) z = - 1 _ 2 ((lnw a ) z + epi (\nw a ) p ) . 
l -+- ep^ i -+■ ep^ 

Let us list some simple solutions. If w a is constant, then also Vi a . If pi = Pi, then 
v. _ h j (pPi)~ 1/2 if w a = p 

Via-*a\ (pp^l + ep-^) if W a = Pa , 
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where p a shall be given by the same expression as some pi, but with in general different constant, 
say a' a . ki a is an arbitrary constant. If pi = pi, then 



Via — kia 



(Pi/p) 1/2 if w a = p 

(Pi/p) 1/2 {1 + epiPa)' 1 if w a =p a 



More complicated solutions are now obtained by noting the following. 

• If Vi a is a solution for w a , then V~} is a solution for w' 1 . 

• If w a is the product of two solutions of (3.11), then V, a is the product of the respective 
solutions for the factors. 

3.3 The Sylvester equation and the solution formula 

Recalling that U, V and Q are ^-independent, the formula for (^-independent) g in (2.8) requires 
X = e 6 ' X with ^-independent X. (2.7) becomes the ^-independent Sylvester equation 

X P -QX = VU . (3.12) 

If spec(P) n spec(Q) = 0, (3.12) has a unique solution, for any choice of the matrices on the right 
hand side. The two matrices P and Q are then independent, hence Theorem 2.1 implies that 8 

g = (I + U(QX)- 1 V)g (3.13) 

solves the non-autonomous chiral model equation (1.1). Obviously, scaling U or V with an arbitrary 
non-zero constant leaves g invariant. We recall from [9] (see Remark 4.4 therein) that 

detP , . „ 

det<7= — det 50 • (3.14) 
det Q 

Example 3.7. Let P be diagonal, as in Example 3.2, and also Q, with eigenvalues given by an 
expression of the same form as pi. If they have no eigenvalue in common, i.e. {p{\ n {<ji} = 0, then 
the unique solution of (3.12) is given by the Cauchy-like matrix 

X ii = ^¥. (3.15) 

Pj ~ Qi 

A vast literature exists on solutions of the Sylvester equation (3.12), more generally with non- 
diagonal matrices P and Q (and not necessarily satisfying the spectrum condition that guarantees 
a unique solution). 

Proposition 3.8 ([32-34]). Let spec(P) n spec(Q) = and 

<£(A) = f> fc A fc 

fc=0 

be the characterictic polynomial of P. Then the unique solution of the Sylvester equation (3.12) is 
given by 

n k—1 

x = -^(Q)- 1 ^<p fe Y^Q k ^ vupl ■ ( 3 - 16 ) 

k=l i=0 



8 We have to choose Q invertible and make sure that the solution X of (3.12) is invertible. See [31,32] for conditions 
that guarantee the latter. 



9 



Remark 3.9. (2.9) takes the form 

X p + p- l X + X z P-Q z X-V z U = 0, 
eX z -(X p - p- l X)P + {Q p + p- l Q)X + V p U = 0. (3.17) 

If we drop the spectrum condition for P and Q, these equations also have to be solved. Otherwise 
they are a consequence of our assumptions (see the proof of Theorem 2.1). In the following, they 
will only be used in the proof of Proposition 4.3 in Appendix B. 

Remark 3.10. Using the above results, we also obtain solutions 4> of (2.2), given by the expression 
in (2.8), which in the case under consideration and via (3.2) takes the form 

e cf>zz + (4> P + p~ x $) P = [4> P + p~ l 4> , 4>z] 

(which corrects a typo in (4.3) of [9]). 
3.4 A reduction condition 

It is of particular interest (see Section 4) to find a convenient condition which guarantees that the 
solution matrix g given by (3.13) is symmetric, i.e. g J = g, where T means matrix transpose. The 
following result is easily verified by a direct computation. 

Lemma 3.11. Let P 1 = -eQ' 1 and g] = g . If V solves (3.10), then U = (V g ) J solves (3.9). 

Proposition 3.12. Let P 1 = -eQ' 1 , spec(P) n spec(Q) = 0, g J = g , and U = (V g ) J • Then 
g given by (3.13) is symmetric. 

Proof. Using g^ = go and U = (V go) J , the Sylvester equation (3.12) and its transpose lead to 

Q [QX - (QXy] = Q [XP - P J X J ] = [QX + eX 1 P~ 1 ] P = [QX - (QX) 1 } P . 

In the last two steps we used P r = —eQ 1 . Now the spectrum condition implies QX = (QX) J . 
Together with g^ = go and U = (V go) 1 , this shows that g given by (3.13) is symmetric. □ 

Remark 3.13. Lemma 3.11 and Proposition 3.12 also hold with transposition replaced by any 
involutory anti-automorphism of the matrix algebra, hence in particular for Hermitian conjugation. 
The Hermitian reduction of the non- autonomous chiral model appears in particular in the context 
of the (electro- vacuum) Einstein-Maxwell equations in four dimensions with two commuting Killing 
vector fields (also see [9]). 

In terms of the matrix 

T := -QX = e{P r )- l X , (3.18) 

which is symmetric under the assumptions of Proposition 3.12 (see the proof of Proposition 3.12), 
the Sylvester equation (3.12) takes the form of a Stein equation, 

T + eP J TP = V g V J . (3.19) 

Implementing the assumptions of Proposition 3.12 in the solution formula (3.13), we have 

g = [I - goV^T- l V) go . (3.20) 
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If P is diagonal, i.e. P = diag(p l5 . . . ,p n ), and if pipj / — e for all then the solution of (3.19) 
is given (via Example 3.7) by 

r (v^k (321) 

This is essentially the corresponding matrix (usually denoted by V) in the Belinski-Zakharov method 
[10,11]. 

Remark 3.14. From Remark 2.4 we deduce the following superposition result. Let (Pi, Vi), i = 
1, . . . , N, be solutions of (3.4) and (3.10) with Q { = -e (Pj)~\ for the same seed solution go- Then 
P = block-diag(Pi,...,Piv) and V = (V\,...,V J N y solve (3.4) and (3.10) with Q = -e^P 7 )^- 
If spec(P) n spec(— e(P T ) _1 ) = 0, and if (3.19) has an invertible solution T, then (3.20) is again a 
symmetric solution of (1.1). 



4 Solutions of the vacuum Einstein equations 

In D dimensions, let us consider a space-time metric of the form 

ds 2 = gij dx l dx^ + / (e dp 2 + dz 2 ) , 

where the (real) components g^, i,j = 1, . . . , m, and the function / only depend on the coordinates 
p and z. The metric then obviously admits m = D — 2 commuting Killing vector fields. 9 For the 
next result, see e.g. [11,14,16,18]. If 

(1) the matrix g = (gij) satisfies det g = —ep 2 , 

(2) g solves the m x m non- autonomous chiral model equation (1.1), 

(3) / is a solution of the compatible system of linear equations 

(In/), = -- + ^tr(U 2 - eV 2 ) , (ln/) 2 = -Ur(WV) , (4.1) 
p 4p 2p 

where U := pg p g~ l and V := pg z g~ X , 

then the metric is Ricci-flat, hence a solution of the vacuum Einstein equations (with vanishing 
cosmological constant). 

Since g has to be real and symmetric, corresponding conditions have to be imposed on the 
matrix data of the class of solutions obtained in Section 3, so that these solutions determine Ricci- 
flat metrics (also see [9]). Such conditions have been found in Section 3.4, and, accordingly, in all 
examples of this section we shall set 

Q = -(p- 1 y, U = (Vg y. (4.2) 

Typically we will regard these equations as defining Q and U in terms of P, V and g$. In most of 
the examples below, P is diagonal. The (symmetric) solution of the non-autonomous chiral model 
is then given by (3.20) with T in (3.21). If P is not diagonal, we have to proceed via the solution 
(3.16) of the Sylvester equation (still assuming that the spectrum condition holds) and (3.18). 

The determinant condition (1) can always be achieved, according to the following remark (also 
see [11,19]). 

9 It is not the most general metric admitting m = D — 2 commuting Killing vector fields. See e.g. [11]. 
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Remark 4.1. By taking the trace of (1.1), one finds that det g and any power of it is a solution 
of (1.1) in the scalar case. We further note that g = w g, with any non-vanishing solution w of the 
scalar equation, is again a solution of (1.1). This shows that, for a given solution g of (1.1), there 
is always a choice of w such that det g = —ep 2 , if m is odd. This is also true if m is even, provided 
that g has the correct signature. 

In the next subsection, we address (4.1). Then we sketch a useful procedure to construct non- 
diagonal metrics from diagonal ones in such a way that the diagonal metric is recovered by setting 
some parameters to zero. A collection of relevant examples in four and five space-time dimensions 
follows. The method is indeed of most interest for D = 4 and D = 5. The higher the number of 
dimensions, the more restrictive is the assumption of D — 2 commuting Killing vector fields for the 
set of solutions of the vacuum Einstein equations. 

4.1 Solutions of the equations for the metric function / 

Example 4.2. For a diagonal solution go of (1.1), so that (go)aa solves (3.11), (4.1) leads to 



fo = Kp 1 Y[ fa , 



a=l 



where k is an arbitrary constant and f a has to be a solution of 

(lnf) p = ^[(lnw) p 2 -e(\nw) z 2 ) , (lnf) 2 = faw) p (\nw) z , (4.3) 

with w replaced by (go)aa- If to is a constant, then also f. Let us write f[w] for the solution of 
the above equations for a given solution w of (3.11), and let oc denote equality up to a non-zero 
constant factor. In particular, we have 

f\p,i] oc ^ where pi = ±y/ [z - dj) 2 + ep 2 - (z - a*) 

(cf. (3.8)). More generally, we find 



--1 j=l i<k 

r s 



(n^'^)nn^^]- 1 , 



j<l i 3 

where p^ is pi with the constant a% replaced by some constant a\ , and we introduced the abbreviation 

PiPj _ Pi Pj 



$[Pi,Pj] 



PiPj + ep 2 2 (at - Qj) ' 



where the last identity holds if o« / aj. The pi (and also the /iQ need not be distinct in the above 
formula. Using these results in our formula for /o, we have the solution of (4.1) for a large class of 
diagonal solutions of (1.1). See Appendix A for some details, and also [10]. 

For the proofs of the following propositions, see Appendix B. 
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Proposition 4.3. Let fo be a solution of (4-1) for a seed go (solution of the Miura equation). 
Let g be a corresponding solution of the non- autonomous chiral model from the family obtained in 
Section 3 (i.e. given by (3.13), where the ingredients are subject to the respective equations). Then 

, _„ det P det Q det X 

f = KfoP T2 T2 , 

[det(J + eP ) det(I + eQ )] x /2 

with an arbitrary constant k, solves (4-1) with the right hand sides evaluated with g. 

Proposition 4.3 does neither assume that go and g are symmetric, nor the reduction conditions 
(4.2). The following corollary now specializes to the case of this reduction. 

Corollary 4.4. Let go be symmetric, Q J = —eP 1 , spec(-P) D spec(Q) = 0, and U = (Vgo) J - 
Then the solution of (4-1) corresponding to g given by (3.20) is 



f = *fop- n (detP) ! 2 detr, 

det(/ + eP ) 



with an arbitrary constant k. 



A convenient formula for the determinant of the solution T of the Stein equation seems not to 
be available in the literature. 

If we have to modify a solution g in order to achieve the determinant condition (1), the following 
result is of great help. 

Proposition 4.5. Let f be a solution of (4-1) for a given solution g from the class obtained in 
Section 3, and let fr w \ be a solution of (4-1) for a scalar solution w of the non- autonomous chiral 
model (1.1). Then 

M^/Hr 1 ^/, (4.4) 

J(detff) 

solves (4-1) with g = wg (which also solves (1.1) according to Remark 4-1)- If det g = —ep 2 , this 
reduces to 

Due to these results, the problem of explicitly computing solutions of the vacuum Einstein 
equations, based on the results in Section 3 and with a diagonal seed, essentially boils down to 
that of solving the linear equations (3.10) for V. For diagonal Q = —eP , solutions of the latter 
have already been obtained in Example 3.6. Although we do not yet have general results in case 
of a non-diagonal Q, corresponding examples are treated below in Examples 4.8 and 4.9. 



4.2 From diagonal to non-diagonal solutions 

In the following, we impose the reduction conditions (4.2), which determine Q in terms of P, and 
U in terms of V and go, and restrict our considerations to the case e = 1. 

We translate a procedure due to Pomeransky [35], which has been formulated and applied 
in the Belinski-Zakharov (BZ) approach, to our framework. We start with a diagonal (hence in 
particular static) solution of the above dimensionally reduced vacuum Einstein equations. The 
corresponding solution g of (1.1) is thus diagonal. According to Remark 2.3 we may still simplify 
it by multiplication with a suitable scalar solution of (1.1). Now some (anti-) solitons are removed 
from it in the following way. 
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• Removal of a soliton p^ = p /x- (i.e. a soliton at z = Oj) from g aa . This means multiplication 
of 9aa by — 2 . In the BZ language, the corresponding trivial BZ vector has a 1 at the ath 
position and otherwise zeros. 

• Removal of an anti-soliton pj = —pj 1 = ~P~ l ^j (i- e - an anti-soliton at z = dj) from gpp. 
This means multiplication of gpp by p~ 2 = —p 2 . 

Afterwards these solitons are reintroduced, with more freedom, in the following way. 

• In order to reintroduce the corresponding soliton, set Pkk = Pi-, for some k. The kth row 
of the matrix V that solves (3.10) (with Q J = —P ), or rather the vector formed by its 
constant coefficients, generalizes the aforementioned trivial BZ vector. It's ath component 
has to be non-zero. 

• To reintroduce the corresponding anti-soliton, set P// = pj, for some I. The Zth. row of the 
matrix V generalizes the corresponding trivial BZ vector if its /3th component is non-zero. 

After multiplication with a suitable scalar solution of (1.1), in order to achieve the determinant 
condition, this leads to a non-diagonal (and thus typically stationary) generalization of the original 
diagonal metric. We recover the latter by suitably fixing the parameters introduced in the second 
step. 

In particular, this offers the possibility to reconstruct a known non-diagonal solution from its 
diagonal specialization (provided the latter can be achieved) via the solution-generating technique. 
Some well-known examples of four-dimensional space-times are recovered in Section 4.3 in this way. 
The above procedure will be applied in the five-dimensional case in Section 4.4. 

Remark 4.6. According to an observation at the end of Remark 2.3, any transformation V i->- 
WV with a constant n x n matrix W that commutes with P J , i.e. [W,-P T ] = 0, leaves g given 
by (3.20) invariant. If P is diagonal, any diagonal W commutes with it. We can then use such 
a transformation of V to scale in each row of V one of the non-zero constant coefficients to a 
fixed value like 1. In the above step of reintroduction of solitons, without restriction of generality 
we can thus fix the constants appearing in the respective positions of the rows of V . If P = P r 
(see Example 3.4), then W is an arbitrary constant, lower-triangular Toeplitz matrix and we can 
again rescale one non-zero coefficient in each row of the corresponding matrix V to a chosen value 
(different from zero), without changing the corresponding solution g. This feature generalizes to a 
P that is block-diagonally composed of matrices P r . 

4.3 Solutions of the vacuum Einstein equations in four dimensions 

Let D = 4, hence m = 2, and e = 1. We thus consider stationary and axially symmetric solutions 
of the vacuum Einstein equations. The metric components <ju, i,j = 1,2, should then refer to 
coordinates x 1 = t (time) and x 2 = ip (angle around the symmetry z-axis), p is the coordinate 
distance from the axis. This interpretation in general imposes additional conditions on the metric. 
The simple solution of (1.1), 

» = ( o 1 2 ) , (4.6) 

corresponds to the four-dimensional Minkowski metric ds 2 = —dt 2 + p 2 dip 2 + dp 2 +dz 2 in cylindrical 
coordinates. The solution /o of (4.1) with g given by (4.6) is simply a constant, which can be set 
to 1 (since this can be achieved by a coordinate transformation). 
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In order to facilitate comparison with the relevant literature, in the following we will mainly 
use the pi introduced in (3.8), instead of the pi, see Notation 3.3. Furthermore, we will frequently 
use the abbreviation (also see [18]) 

fcij = p 2 + m Hj = 2 (a, - a,j) ^ j , (4.7) 

Hi jij 

where the last identity holds if a« / aj . 

Example 4.7 (Kerr-NUT). In order to recover the Kerr-NUT metric, we start with its static 
specialization, the Schwarzschild metric. The latter corresponds to the following solution of (1.1), 



ill o 
p 



~9=\ n Ml „2 Mi ) ' (4-8) 



/'2 



Next we remove a soliton at z = a\ from gu and a soliton at z = ai from $22- The resulting matrix 
is simplified by rescaling it with the inverse of 

w = p~ 2 P1P2 ■ 

We obtain the seed solution 

go = w~ 1 gdi&g(-(p/p 1 y 2 ,-(p/p 2 y 2 ) , 
which is nothing but (4.6). Next we reintroduce the two solitons (n = 2) via 

P = diag(A^). 
The respective solution of (3.10) (with Q = —P 1 ) is 



V 



Vll Vup 1 1 
V21 V 2 2 P2 1 



with constants Vi a (see Example 3.6). According to Remark 4.6, we can set v\\ = V22 = 1 without 
restriction of generality. Let us rename the remaining parameters in V, 

ci = V12 , c 2 = v 2 i ■ 

Now we obtain the matrix T from (3.21). Then g = wg, with g given by (3.20), satisfies the 
determinant condition, i.e. det 5 = — p 2 , and has the components 

5ii = ^ ((c? + c 2 ) p 2 pip 2 (pi ~ P2) 2 ~ P\P2 ((1 + c\c%) 1Z\ 2 - 2ci 02^11^22) ) , 
912 = (Ml ~ ^ Ul2 (ci (c 2 p 2 - p 2 ) K U + c 2 (c 2 p 2 - p 2 ) K22) , 

522 = \ ( - ^ = 112)2 (cj p 8 + c 2 p\pf) + ((pi + c 2 c 2 pj) n 2 2 - 2 ClC2 p\p\ n u n 2 2) ) , 



where 



h = (pi - p2) 2 (c\p A + clplpl) + (pi + c\c\ p 2 2 ) Tl\ 2 ~ 2c i c 2 P1P2 Tl\{R22 ■ 
Furthermore, according to Example 4.2, we have 

P\P2 



f(w) 



oc 



\ITl\\R22 Tl\2 
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Corollary 4.4 and Proposition 4.5 then yield 



with a constant k. Setting a 2 = —a\ = cr, 

ci = m + r °" , c 2 = fl + ^ , k = C — p- , a = \/m 2 - o 2 + b 2 , 
a + b m — a 4a z 

with new constants a, b, m and C, and passing over from the coordinates p and z to new coordinates 
r and 6* via 



p = y/ (r — m) 2 — a 2 s'm9 , z = (r — m) cos 6* , 

we obtain the Kerr-NUT metric in Boyer-Lindquist coordinates, as given by (8.48) and (8.49) in 
[11]. We can compose N Kerr-NUT data blockwise into larger matrices (the new P and V are 
then 2N x 2N, respectively 2N x 2 matrices), and obtain again a symmetric solution of the non- 
autonomous chiral model equation, see Remark 3.14. The determinant condition is again achieved 
by multiplication with a suitable scalar. In this way we can recover the multi-Kerr-NUT solutions 
[36]. 

Choosing for P the non-diagonal solution of (3.5) resulting from an n x n Jordan block matrix 
A n (see Example 3.4), we should expect that the resulting family of solutions of the non-autonomous 
chiral model can be obtained alternatively via a "soliton coincidence" limit ("pole fusion" in the 
Belinski-Zakharov formalism [11]) of the family of solutions obtained with a diagonal P (with 
distinct eigenvalues). This is confirmed by the following example. 

Example 4.8. Again, let n = 2 and go as in (4.6). Now we choose 

p ( 1 



P = 

Mi 



1 



where fii = 9\ — z with 9! = ^Jz 2 + p 2 . This is P 2 in Example 3.4 (with a = 0). We find the 
solution 



V 



vu v 12 p 1 1 

V21 V 2 2 Pi 1 ~ 2vu TZii 



of (3.10), with constants V{ a . Computing the solution of the Sylvester equation via (3.16), the 
symmetric matrix T is obtained from (3.18). We have to rescale the resulting g, given by (3.20), to 
g = w g with w = p~ 2 pf, in order to arrange the determinant condition. It turns out that g only 
depends on the parameters Vi a via a = vu/vn and f3 = det(t>j a ) / 'v\ x . Furthermore, we find that g 
also results from the corresponding (Kerr-NUT) solution in Example 4.7 in the limit a — > 0, after 
setting c\ = a and c 2 = a~ 2 (a — 2(3 a). 

The Kerr-NUT metric consists of two solitons of multiplicity one. The Tomimatsu-Sato metrics 
[37] are known to generalize it to two solitons of multiplicity S > 1 (see e.g. [II]). 10 We should 
then expect that, for fixed S, this solution can be obtained alternatively with P consisting of two 
blocks Ps as given in Example 3.4 (corresponding to the 6 x S Jordan block As). For 5 = 2 this is 
confirmed in the next example. 



The Tomimatsu-Sato metrics are in fact limiting cases of multiple Kerr-NUT solutions [36,38]. 
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Example 4.9 (Tomimatsu-Sato). In the static limit, the 5 = 2 Tomimatsu-Sato metric (also see 
[39]) reduces to the Zipoy-Voorhees metric. The latter corresponds to the diagonal solution 

g = diag (-^ 2 p 2 2 , p 2 p\p 2 2 ) 

of (1.1). Now we remove twice a soliton at z = a\ = a from gu, and also twice a soliton at 
z = 0,2 = — cr from #22- With a simplifying scaling we obtain the seed 

— 1 ~ t I —4 4 —4 4\ —4 2 2 

50 = w ffdiag(p ^ 1; p p, 2 ) , w = p p x p 2 . 
The resulting go is again (4.6). Next we reintroduce the two double-solitons via 11 



/ 


PP1 1 


-2 P n^ 















PPi 1 


















PP2 1 


-2 P n 22 l 




V 











PP2 1 


) 



The solution of (3.10) is then given by 



V = 



\ 



V21 v 2 2 Pi 1 - 2v 12 Till 

\ f 4 l V42 P-2 1 ~ 2v 32 1^22 I 



Again, the matrix T is obtained via (3.16) and (3.18). Next we have to compute g = wg, with g 
given by (3.20), and 



-6 --2 (detP) 2 detr 3 
f = Kp wf {w) — _^ where /(«,)= p 



det(I + P 2 ) 



2 2 



^Il7£l2^22 

According to Remark 4.6, without restriction of generality we can set 

vu = V21 = V32 = V42 = 1 • 

Choosing 



V12 = V22 = V31 = V41 



1 + P 



k = a 



1 + P 
4p 



in terms of prolate spheroidal coordinates x, y, given by 

p = ay/(x 2 _ i)(i_ y 2) j z = axy ^ 



(4.9) 



and after a coordinate transformation £ i-> i — 4<rp 1 q( / 9, we obtain the (5 = 2 Tomimatsu-Sato 
metric 



ds 2 = --dt 2 + 8a(l-y 2 )?—dtdcp + a 2 (l-y 2 )-^-d(p 2 + 



D 



B 



^ + ^ 



p 4 (x 2 - y 2 ) 3 \x 2 - 1 1 - y 2 



1 Noting that pi/SHi = 2p/1Zn, this consists of two blocks of the form P2 in Example 3.4. 



17 



with 



A = (l-y 2 ) 4 g(Z) 9 (-Z), 

B = ( (x) + qV) 2 + 4qV (p x 3 + 1 - (p x + 1) y 2 ) 2 , 

C = q 2 (px + l)/(-y 2 + 3) + (-2q 2 (p x 3 + 1) + (p x + 1) g(x)) y 2 - (2p x 3 - p x + 1) g(x) , 
D = p 2 q 4 (x 2 -l)y 8 + 4q 2 (px + l)(p 3 x 3 + 3p 2 x 2 -p 3 x + 4px-3p 2 + 4)y 6 

-2q 2 (3p 4 x 6 + 4p 3 x 5 - 3p 4 x 4 + 8p 3 x 3 + 37p 2 x 2 - 12p 3 x + 48p x - 13p 2 + 24) y 4 

+4q 2 (p 4 x 8 - p 4 x 6 + 4p 3 x 5 - 4p 3 x 3 + 15p 2 x 2 + 24p x - 3p 2 + 12) y 2 

+g(x) (p 4 x 6 + 6p 3 x 5 - p 4 x 4 + 16p 2 x 4 - 12p 3 x 3 + 32p x 3 + 15p 2 x 2 + 6p 3 x - 15p 2 + 16) , 

where p 2 + q 2 = 1 and 

g(x) = p 2 x 4 + 2px 3 -2px-l, Z 2 = * - \ . 

1 - r 

These are the expressions (32)- (35) in [39]. 12 

Remark 4.10. The m = 2 non-autonomous chiral model is also related to the D = 4 vacuum 
Einstein equations with two commuting Killing vector fields in another way, via the Ernst equation 
(see e.g. [12]). In this case, one has to impose a different determinant condition: detg = 1. 
This allows a constant seed go, and thus constant matrices U, V, in which case the application of 
Theorem 2.1 is considerably simplified. The solutions obtained in this way include the multi-Kerr- 
NUT metrics (see [9] for details). 13 A similar construction with m = 3 leads to the (electrically 
and magnetically) charged generalizations, the multi-Demianski-Newman metrics [9]. 

4.4 Solutions of the vacuum Einstein equations in five dimensions 

Let m = 3 and e = 1. In this case it seems that no relevant solutions can be obtained by choosing as 
the seed a diagonal solution of the non-autonomous chiral model corresponding to five-dimensional 
Minkowski space-time. It then becomes a subtle problem to choose a suitable seed solution. Here 
insights about the "rod structure" of the putative axis p = are of great help [11,15,35,42,43]. 
Some important solutions that have been obtained or recovered previously in the Belinski-Zakharov 
approach will now be presented in our framework. We use the procedure outlined in Section 4.2. 

4.4.1 Myers-Perry black holes 

The higher-dimensional generalization of a static black hole is given by the Schwarzschild-Tangherlini 
solution [44]. The corresponding diagonal solution of (1.1) is 

g = diag( ,fi2, — ) • 

M2 Ml 

The seed go is obtained from g as follows: 

• Remove a soliton at z = a\ from 533. 

• Remove an anti-soliton at z = a<i from 522- 

12 We obtained a factor q 4 in the expression for D, whereas there is a factor q 2 in [39]. We believe that the latter 
is a typo. 

13 In [40,41], solutions of the five-dimensional Einstein-Maxwell equations are constructed from a pair of solutions 
of the m — 2 non- autonomous chiral model with symmetric g and det g = 1. 
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• Multiply the resulting matrix by w 1 , with w = — /xi/// 2 , to achieve a simpler form. 
This results in 

P 2 

g = diag(l, — ,// 2 ) • 
Mi 

According to Example 4.2, the associated solution of (4.1) is (up to a constant factor) 

, M2 
Jo 



Next we reintroduce the soliton and the anti-soliton via 

P = diag(A-^). 
Ml P 

The first row of the 2x3 matrix V should have a non-zero entry in the third component, which 
means we recreate the soliton at z = a\ in the respective diagonal component of the seed. The 
second row of V should have a non-zero entry in the second component, which means we recreate 
the anti-soliton at z = a 2 . According to Example 3.6, the solution of (3.10) (with Q = —P ) is 
given by 

V = 



V21 V 22 V 2 3 ^ 



with constants Vi a . The corresponding solution of the Stein equation (3.19) is obtained from (3.21). 
With g given by (3.20), the new solution 

g = wg 

satisfies the determinant condition: det g = —p 2 . We find 

n 12 



f(w) 



p^TZ n TZ22 



(up to a constant factor). Without restriction of generality, we can set v\ 3 = V22 = 1 (see Re- 
mark 4.6). With the restrictions and redefinitions 

V12 = V 2 3 = 0, Ci = Vu , C 2 = V 2 1 , 

use of Corollary 4.4 and Proposition 4.5 yields 

f = Kfl2 c i c 2 Mi M2 - (Mi ~ M2) 2 (mi (C1M1M2 + 4p 2 ) + 1Z 2 12 ) 

(Mi - M2) 2 ^11^-12^-22 

with a constant k. The corresponding metric is given by 

ds 2 = 511 dt 2 + g 2 2 dip 2 + 533 dip 2 + 2 512 dtdip + 2 g 13 dtdip + 2 g 23 dip dip + f (dp 2 + dz 2 ) . 
Setting 

a± = —a , a 2 = , 
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with a constant a, in terms of prolate spheroidal coordinates x,y, given by (4.9), we have 

g u = - (4a x + (aj - a 2 2 ) y - p 2 ) u- 1 , 

912 = -aiPo (l - y)^ 1 , 

gn = -a 2 po (l + y) u' 1 , 

922 = (4ax + a 2 -a 2 + p 2 + 2a 2 p 2 (l-y)u J - 1 ) , 

923 = ^aia 2 pl (1 - y 2 )^" 1 , 

533 = (4ax-a 2 + a 2 + p 2 + 2a 2 p 2 (l + y)uj- 1 ) , 

f = ^, u = lax+{a\-al)y + pl. 



Performing a linear transformation of the coordinates t, if, ip, 



i \ / . , ai a 2 c 2 oi ci a 2 

{t, <p, ip) >->• {t - oi ip - a 2 ip , ¥> + - — ip , — — <f ip), 

c 2 4 a 4cr ci 



and setting 



4a (pg - o? + a| - 4a) 2 _ 4a (p§ + a? - a 2 . - 4a) 

c 2 



PO + fl l " a 2 + Po " a l + a 2 + 4CT 



we recover the analogue of the Kerr metric, i.e. the Myers-Perry solution [15, 45, 46] of the five- 
dimensional vacuum Einstein equations, precisely in the form of equation (19) in [35]. Switching on 
the constants vu and V23 (which we set to zero above), leads to a more general class of space-times. 

4.4.2 Black saturn 

A black saturn [47] (also see [48-50]) is a black hole surrounded by a black ring. 14 In the following, 
we show how the black saturn solution, originally obtained in [47], can be recovered in our approach. 
Let us start with the same static solution of the non-autonomous chiral model as in [47], 

~ , MiP3 P 2 Pi P2P5, 
g = diag( , , ) . 

P2P4 P3P5 Pi 

This is motivated by a rod structure analysis (also see [15,42,43]). A suitable seed solution is then 
constructed as follows. 

• Remove anti-solitons at z = a\ and at z = 03 from gu. 

• Remove a soliton at z = a 2 from gu. 

• Multiply the resulting matrix by w~ x with w = to simplify its form. 



This results in 



-1 ~v 1 P f4 1 -n i - / PlP4 P3P5 , 
go = w #diag( 1, 1) = diag(— , , — =-) . 

PlPl P4 P2P5 P 



14 Black rings are similar to black holes, but with horizon topology S 1 x S D 3 . They only appear in D > 4 
dimensions. See [51-69]. 
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The corresponding solution of (4.1) is obtained via Example 4.2, 

,. , 2 ^-12^-15^-24^-45 
JO = « M4 A*5 



^-22^-33 ^-14^-25^-35'^-44'^-55 ' 

with a constant fc. We have n = 3 and choose P diagonal with 

P\l = -Hl/P, P22 = P/^2, P33 = -M3/P- 

The solution of the equations for V is then given by 

„, „. M2M5 TlwTili „, Ml gggj^gg 

^21^7 ^22 w/mW22 7J25 ^23 M ^ 3M5 

\ U 31 M3 ^32 TC 13 ^3 4 ^33 ^33^35 

with constants Vi a . Without restriction of generality, we can set 

vu = V21 = V31 = 1 . 

The case considered in Section 2.2 of [47] should then correspond to the subclass of solutions given 
by 

V12 = V22 = v 33 = . 
The authors of [47] then only elaborate the special case 

V32 = 

further. In this case, however, the removal of the anti-soliton at z = 03 from the original static 
metric and the subsequent reintroduction via P33 is actually redundant. This means that the 
black saturn space-time can already be obtained from n = 2 data. The solution T of the Stein 
equation (3.19) is given by (3.21). With g given by (3.20), we set g = w g to satisfy the determinant 
condition. Then g reduces to g if v 13 = v 23 = 0, as expected. We find 

^-12^-23 

J(w) 



V / ^-ll^-22^-33 7^13 

(up to a constant factor) and, using Corollary 4.4 and Proposition 4.5, 

i -6 - ,- 3 (detP) 2 detr 

f = -Kp Wf f (] , 

( ' det(/ + P ) 
which turns out to be a lengthy expression. Setting 

(oti - a2) 2 (a 2 - a 3 ) 2 



v 13 = 2ci (ai - a 5 ) , t> 2 3 = c 2 [4(a 2 - a 3 )(a 2 - a 4 )] 1 , k = 4- 



(ai - a 3 ) 2 



with constants Cj, we obtain the metric 

ds * = -^ dt+ ^- di; y + H l (k 2 v {dp 2 + dz 2 ) + g v + g # 2 ) , 
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where 



= (M + cfMi + C2M 2 + C1C2 M 3 + c 2 c| M 4 ) , 

H 2 = ^(^M -cl^M 1 -cl^M 2 + c 1 c 2 M 3 + clcl^M l 



F = 



5 

PlPb{Pl - /"3) 2 (/"2 - /U4) 2 (fJ^jj) ^13^-14^23^-24^-25^-35 



i=l 



with 



M) = ^2/^5 (/"I - M3) 2 (>2 - M4) 2 T^\2^A^h , 

Mi = p 2 p\p 2 P3PAPb (Pi ~ P2) 2 {p\ ~ Pb?{P2 ~ PA? R-23 > 

M 2 = p 2 plp 2 ~ 1 P3P4P5 (Pi - P2f(pi ~ P3) 2 ^-14^-25 > 

M 3 = 2 / U 2 /i 3 /i 4 /i 5 (jLti - /X 3 )(//i - /U 5 )(^2 - /U4) ^-11^-22^-14^-23^-25 , 

M 4 = /if/I^plpl (Pi ~ P5f ^- 2 2^-25 > 

2 = (ciJRi ^M M 1 - c 2 y\ 2y / M M 2 + c 2 c 2 9W^i^4 - cic^i V-M2M4) , 



Co' 



and 9tj = y/(z — a^) 2 + p 2 . This is the black saturn metric 15 as given in [18,47], with some obvious 
changes in notation, but some deviations in the factors \ii in the expressions for M 2 , M3 and M4. 

4.4.3 Double Myers-Perry black hole solution 

In order to recover the double Myers-Perry black hole solution obtained in [70], we start with the 
matrix that determines a static two-black hole solution (cf. (3.1) in [71]), 

_ ,. , PlPA P 2 P3 P2P5, 

g = diag( , , ) . 

M2A*5 PlPA P3 

Removal of two solitons p 2 ,p5 and two anti-solitons P\,Pa from gu, and simplification with a 
suitable factor, leads to 

-1 - A- ( ( Plpl ^ 1 1 A a- ( 1 P 2 P3 P±P±\ P2P5 

g = w flidiagl - (-2-3J , 1, 1 = diag(-l, , ), w = . 

v K P 2 PV 1 P2Pb P3 PlPA 

According to Example 4.2, the corresponding solution of (4.1) is given by 

f Q = Z 21 ^ 4 ^ TJ t/r~^ ~* ^13^23^-34^-35 



P3 ^ 7 = j ' ^-14^-25^-33 

We have n = 4 and reintroduce the removed solitons and anti-solitons via 

P = diag(-p~Vi, - P' 1 Pa, p P 2 l , p P5 1 ) ■ 



15 Actually, additional conditions have to be imposed on the remaining parameters in order to achieve asymptotic 
flatness and absence of naked and conical singularities, see [47,49]. 
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The solution of (3.10) is given by 



V21 

V31 V32 



^12^15 
v 7^24^45 

U22 P 2 n 3 



•34 

M3 ^22^25 
UoM r > 7?-3K 

v U41 V42 At,/;';, 



u i3 ^i^n 

„, M4 73-34 

"23 7e 14 7e 44 



W33 



«43 



M3 73.127^24 
M1M2M4 7^23 



M3 Tjjjj'gjs 



JU1//4A15 



^35 / 



In order to recover the double Myers-Perry black hole solution presented in [70], we reduce the set 
of solutions obtained in this way by restricting V to 



V 



/ 1 
1 
1 

V 1 








\ 


J 



We immediately notice that this means in particular "trivializing" the solitons at z = a>2 and z = 05 . 
Hence, the solution obtained with this special choice of V can already be obtained from n = 2 
data. We should expect, however, a 4-soliton transformation to be necessary in order to generate a 
(sufficiently general) double black hole solution, which suggests to explore the above more general 
solution. This will not be done here, and we return to the special case with the above restricted V. 
Again, we obtain T from (3.21). Let g be the resulting solution (3.20), multiplied by w to achieve 
the determinant condition. From Example 4.2, we obtain 







i^3 



1 7^12^.15^-24^-45 
^-14^-25 



and then / via Corollary 4.4 and Proposition 4.5. Setting 

b(ai - a 3 ) 



V\2 



2(ai - a 2 )(ai - a 5 ) ' V22 2(a 2 - a 4 )(a 4 - a 5 ) ' 
(ai - a 2 )(ai - a 5 )(a 2 - a 4 )(a 4 - a 5 )\2 



c (a 3 - a 4 ) 



i?i V H 2 / P2P5 H 2 P3 F 



(ai - a 4 )(a 2 - a 5 ) 
with constants 6, c, this results in the metric 

P2P5 ^2 A*3 



where 

Hi 
H 2 



M + 6 2 Mi + c 2 M 2 + 6cM 3 + 6 2 c 2 M 4 , 
P 2 (PlPi 



P2P5 





f = Mi(^-^4) 2 (n^)- 



M4 

^12^14^15^24^25^45 



M - 6 2 ^Mi - c 2 ^M 2 - fecM 3 + 6 2 c 2 ^- M 4 ) , 
pi Pi p\P4 > 



i=i 



^13^-23^-34^35 



= -2(-^-) (fe£Hi VM Mi + c^VMqMs - 6 2 c^ 4 VMiM 4 - 6c 2 <Ki y / M 2 M^j , 



<P2P5 
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with 



m = M2M3M5 (mi - M4) 2 nj 2 nj 5 nl 4 nl 5 , 

M l = M1M2M3M5 (Ml - /"3) 2 ^14^24^45 > 

m 2 = M2M3M4M 2 (M3 - M4) 2 nl 2 nl 4 nl 5 , 

M 3 = 2^i/i^/i 3 /i 4 ^^ (/Xi - jU 3 )(/X 3 - /X 4 ) ^11^12^15^24^44^-45 , 
M 4 = P 4 l4vll4vl (Ml - M3) 2 (Ml - M4) 2 (M3 - M4) 2 • 

With obvious changes in notation, this is the metric obtained in [70]. 

Remark 4.11. It is plausible that one can start with the diagonal solution of the non-autonomous 
chiral model corresponding to a static triple black hole space-time (see (4.1) in [71]) and construct a 
space-time with three Myers-Perry black holes. This procedure should continue to produce solutions 
with an arbitrary number of rotating black holes. 

4.4.4 Bicycling black rings 

Let us start with the solution 

- ■,. / MlM5 P 2 M3M7 M2M4M6^ 
9 = diag( , , ) 

M3M7 M2M4M6 MlM5 

of (1.1), which corresponds to a static metric. Removal of a soliton at z = a-j and an anti-soliton 
at z = a\ from gu, and a rescaling, leads to the seed metric 

/M? 1 u a- /M5 M 2 MlM3 M2M4M6x M7 

5o = ^ 5 diag(-^,l,l) = diag(— , , ), w = , 

Mi M3 M2M4M6 M5M7 Mi 

with the following solution of (4.1) (up to a constant factor), 

j. _ M2M4M6 ^-12^-14^-16^-23^-25'^-27^-34 ^-35^-36^-45 ^-47^-56^-67 
M5M7 V^n'^-13^-22'^-24^26^33'^-44'^-46^-55'^-57'^-66 V^-77 

We have n = 2 and shall set 

P = diag(-p"Vi>MM7 1 ) • 
The solution of (3.10) is then given by (see Example 3.6) 



V = 



Vn Bs. Vl2 Kif 14 K 16 gKgRg 

±1 K-15 /il/-Cll/<.l3 ±0 I<12 /-C147C16 

v M3^-57 ^2^4^6^17^-37 Ms'^-27^47'^-67 

21 ^5^37 22 Ml^3^-27^47'??-67 23 M2M4M6 ^-57^-77 



Without restriction of generality, we can set uu = W21 = 1. But we do restrict the class of solutions 
by setting 

V12 = v 2 3 = . 

Again, the solution T of the Stein equation is obtained from (3.21). The resulting solution (3.20) 
of the non-autonomous chiral model has to be modified to g = w g, with w as given above, in order 
to achieve the determinant condition. We have (disregarding a constant factor) 

, _ _i T^n 

1{w) ~ 9 Vk^W 7 ' 
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and, from Corollary 4.4 and Proposition 4.5, 

2 -5 , ,- 3 (detP) 2 detr 

f = K p Wfof ( ) , 

( ' det(J + P ) 
which results in a lengthy expression. Setting 

(a 7 - a 4 )(a 7 - a 5 )(a 7 - a 6 ) . , 2 

«13 = Ci , V 2 2 = 2 -, T7 , K = 4 (a 7 - fill) , 

(a 7 - ai)(a 7 - a 3 ) 2 

we obtain the metric 16 

ds 2 = - _2 f d$ _ dy, _ dij,) 2 + -L ( Gi dp 2 + G 2 # 2 - 2 J (d/> 2 + (fc 2 ) 

-rc 1 V -n 2 -n 2 / il 2 ^ ' 



where 



H-i = M + c 2 Mi + 6|M 2 - &2C 2 M 3 



o- P5/P1 2 P A/f ,2PlP7 ,2 2P7 ,, \ 

iJ 2 = — — M — Mi - 5 2 — 5— M 2 - ft 2 c i — M 3 , 

P3 V P7 PlP7 P 2 Pi ' 

P2P4P6 V A*i Hi ' 

G 2 = ^^(Mo + clM^b^Mi + btc^Ms), 

P3P7 V P 2 P l J 

-i = ^2^f^^V /2 (yM^-c 2 ^yMM3), 

P7 VP2P4P6/ ^ Pl ' 

Pl VP P3P7/ V P / 

J = 6 2 Clp 2 / Ui / U 2 P3P4pip6 (P3 - P7) 2 (P4 ~ P7)(P5 ~ P7)(P6 ~ X 

TZiiTZi2lZi^TZ\4^{IZi^lZijTZ27'Tl7j , 



and 



M = ^4pip6P7(p3-p7) 4 ^ 2 2^ 2 3^ 2 4^16^17^27 ) 

M l = P 2 PlP2P3P 4 P5P6 (Pl - P7) 2 (P3 - P7) 4 ^-15 ^17^27 > 

M 2 = /) 4 /ilP2P3P5P7(P4-Ai7) 2 (P5-P7) 2 (P6-p7) 2 '^ 2 2'^13'^14'^16i 

^3 = P 4 PiP2P3P4P6(P4-P7) 2 (P5-P7) 2 (P6-P7) 2 ^- 4 5^17' 

p P2 ^-23^-25^-34 ^-35^-36 ^-45 ^-47^-56^-57^67 

MlPsP7(P3 " P7) 4 ^12^13^14^15^16^17^24^26^27^37^46 ' 

With obvious changes in notation, this is the "bicycling" black bi-ring solution obtained in [63] 
(also see [64]), except for the fact that we have a minus sign instead of a plus in the expressions 
for u>i and uj 2 - 



5 Final remarks 

We presented a general formulation of binary Darboux transformations in the bidifferential calculus 
framework. Whenever a PDDE can be cast into the form (2.2) or (2.3), Theorem 2.1 can be applied 

16 Here we used a, — ctj = (fn — pij)(p 2 + Him) /{2jniij) to eliminate en — cij. 
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and it will typically generate a large class of exact solutions. Meanwhile a bidifferential calculus 
formulation is available for quite a number of integrable PDDEs. 

We elaborated this general result for the case of the non-autonomous chiral model, considerably 
extending previous results in [9]. We also presented conditions that, imposed on the matrix data 
that determine the general class of solutions, guarantee that the resulting solution of the non- 
autonomous chiral model is symmetric (or Hermitian). If the solution is also real, then it is known 
to determine a Ricci-nat metric, i.e. a solution of the vacuum Einstein equations, dimensionally 
reduced to two dimensions. We essentially solved the equations resulting from the assumptions in 
Theorem 2.1 in the case of a diagonal seed metric, though not yet the V-equations in sufficient 
generality if P is non-diagonal (but see Examples 4.8 and 4.9). All this provides a working recipe 
to compute quite easily solutions of the vacuum Einstein equations. In particular, in the four- 
dimensional case we recovered (multi-) Kerr-NUT (in a different way than in [9]) and the 5 = 2 
Tomimatsu-Sato solution. In the five-dimensional case we recovered single and double Myers- 
Perry black holes, the "black saturn" and the "bicycling black ring" solutions. The more general 
solutions still have to be explored. In view of the complexity of the latter solutions, it is certainly 
an advantage to have now an independent method at our disposal to derive, verify or generalize 
them. Surely further important solutions of Einstein's equations in D > 4 space-time dimensions 
can be recovered using this method and there is a chance to discover interesting new solutions. We 
concentrated on examples in the stationary case e = 1, but developed the formalism as well for the 
wave case e = — 1 [11, 19]. It is not difficult to recover relevant examples in this case too. 

The recipe to construct solutions of the non-autonomous chiral model and the dimensionally 
reduced vacuum Einstein equations, obtained from Theorem 2.1, is - not surprisingly - a variant, 
a sort of matrix version, of the well-known method of Belinski and Zakharov [10, 11, 19]. 

One should look for suitable ways to spot physically relevant solutions within the plethora of 
solutions we now have at our disposal. How are desired properties of solutions, like asymptotic 
flatness, absence of naked singularities and proper axis conditions encoded in the (matrix) data that 
determine a solution? Here the rod structure analysis [11,15,42,43], developed for the Belinski- 
Zakharov approach and frequently used, is of great help. 

Section 3 also paved the way toward a treatment of other reductions of the non-autonomous 
chiral model, which e.g. are relevant in the Einstein-Maxwell case and supergravity theories. 

In this work we only elaborated Theorem 2.1 for a particular example of an integrable equation 
in the bidifferential calculus framework. Although we already applied a more restricted version 
of it previously to several other integrable equations, it will be worth to reconsider them and also 
further equations with the much more general solution-generating tool we now have at our disposal. 

Appendix A: Addendum to Example 4.2 

From (4.3) we find that 

f[p] oc p 1 / 4 , f[w k ] « f [wf , fKw 2 ] (x f[ Wl ] f[w 2 ] d[wi,w 2 ] , 
where $\wi,w 2 ] has to solve 

(\n$[w 1 ,w 2 ]) p = !((lnu/i) p (In u> 2 ) p - e (lniui) 2 (lnw 2 ) 2 ) , 
(\n$[w 1 ,w 2 ]) z = ^((lnu>i) p (lnw 2 ) z + (lnu/i) 2 (lniu 2 )p) 
(also see [64]). It is easy to verify that 

r s 

$[wi ■■■w r ,w' 1 --- w' s ] x Yl n $[Wi,Wj] ■ 

i=l j=l 
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In particular, ^fu^u^] oc *$[wi, W2\ kl ■ Furthermore, we have 

$[p k w u p l w 2 ] oc p kl / 2 w[ /2 w k 2 /2 3[w 1 ,w 2 } . 

It follows that 

f [p k W] OC p k ' 2 / 4 W k / 2 f[w} , f [Al ■ ■ ■ At] OC ( f[ f[pk]) ( II » fa]) > 

fe = l i<jr 

and 

fr Al'-'Atr , f[Al • • • At] flftl • • • Mr] 

from which the main result in Example 4.2 is easily deduced. 

Appendix B: Some proofs 

Proof of Proposition 4.3. Using (lndetY)p = ti{Y p Y~ l ) for an invertible (and differentiable) 
matrix function Y, we obtain 

(In/), = (In J b) p + tr(- p^I + X p X~ l + P p P~ l +Q p Q 1 - eP p P{\ + eP 2 )' 1 
-eQpQil + eQ 2 )- 1 ) , 

(ln/) 2 = (ln/ ) 2 + tr(x,X _1 + P z P~ l + Q Z Q 1 - eP z P(l + eP 2 )~ 1 - eQ z Q(l + eQ 2 )- 1 ) , 

for the expression of / in Proposition 4.3. In order to verify that (4.1) holds, we have to show that 
these expressions equal the corresponding right hand sides of (4.1), evaluated with 

u = -p4> z = -p[4>o,z + (ux~ 1 v) z ] 

= U -p (u z x^v - ux~ l x z 5C x v + UX~Wz) , 

V = V + e[p{U p X~ l V -UX^XpX^V + UX~ 1 V P ) + UX~ 1 V], 

where Uq = —p(j)o,z and Vo = e(pcbo,p + 4>o)- 17 We work on the right hand sides of (4.1) and, as 
intermediate steps, we consider tr(£Y 2 — Uq), tr(V 2 — Vg) and tr(UV — UqVo) separately. First we 
eliminate derivatives of 4>o with the help of (3.9) and (3.10). Then we eliminate VU using the 
Sylvester equation, and V z U, V p U via (3.17). Rewriting (3.17) with the help of (3.12), we obtain 
a version that allows us to also replace all occurencies of VU Z and VU p . Several times one has to 
exploit the cyclicity of the trace. Finally we use 

P p = p- l P(I-eP 2 )(I + eP 2 )-\ P z =2ep- l P 2 {I + eP 2 )- 1 

(which follows from (3.4)), and the corresponding equations for Q, to show that (4.1) holds. 



This is quite a tour de force and a more elegant proof would be desirable. 
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Proof of Proposition 4.5. Using U = U + p (In w) p I and V = V + p (In w) z /, we have 

ti(U 2 -U 2 ) = 2p(\nw) p txU + mp 2 (\nw) 2 = 2p 2 (lnw) p (lndet#) p + mp 2 (\nw) 2 , 
and a corresponding expression for tr(V 2 — V 2 ). Furthermore, 

tv(uV-UV^ = p 2 ((In w) p {lndetg) z + (lnw) z {lndet g) p + m (In w) p (lnw) z ) . 

With the help of (4.1) we obtain 

P P / 2 2\ 

= 2 ^ ( ln det 9 ^p ~ e ( ln w ^ z ( ln det ^ + m 4 i^ ln ^ ~ e ( ln ^ 2 J ' 



1 

fJz " V 



= ^((lnw) p (lndetg) z + (In w) z (ln det g) p ) + m^(lnw) p (lnw) z . 

Next we use 

P ( 2 2\ P 

(Mpf(w)))p = ^[(^w) p -e(lnw) z J, (ln(pf (w) )) z = -(lnw) p (lnw) z , 
and the corresponding equations for /(det 9) an d f\ w det 3)) to deduce that 

(in { {WA T 9) ) = ?-((\nw) p (\ndetg) p -e(\nw) z (\ndetg) z ) , 

V P /(w) /(det g) ' P ^ 

(in { {wd f 9) ) = ?-((\nw) p (\ndetg) z + (\nw) z (\ndetg) p ) . 

V P J(w) /(det 9 ) ' z 2 

Inserting the last expressions in our previous results, we obtain (4.4) by integration. Let us now 
assume that det g = —ep 2 , and thus det g = —ep 2 w~ m . With 

f {p 2 wk) oc w k (pf {w) ) k2 

(see Example 4.2), (4.5) results from (4.4). 
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